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Oscillations of local density of states generated by a single scalar impurity potential are calculated 
for one-dimensional systems with dynamically generated charge or spin gap. At zero temperature 
the oscillations develop at finite wave vector (tt for the Mott insulator and 2kF for ICDW/SC) and 
at frequencies larger than the soliton spectral gap m. Their amplitude has a broad maximum at 
uj ~ 3m, where m is the gap magnitude. 



Modern Single Electron Tunneling Microscopy (STM) 
techniques provide a plephora of data by measuring 
minute details of coordinate dependence of the differen- 
tial conductance ct(x) = dl /AV . An especially rich pic- 
ture emerges in the presence of impurities which create 
Friedel oscillations in cr(x). An ultimate goal for a theory 
is to decipher the interference pattern of these oscillations 
and extract from it information about the spectrum and 
interactions in the underlying system. Since the conduc- 
tance is directly related to the single electron density of 
states (DOE), theory should aim at calculation of this 
quantity: 

a(V^,x) =C(x)p(a; = ey,x); 

p(w, x) = -'imG^^^ [lo] X, x) (1) 
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where C(x) is an a priori unknown matrix element. 

In majority of systems studied by now by STM the dis- 
order is sufficiently weak for momentum to remain a good 
quantum number. For instance, observing sharp peaks 
in the Fourier images of the STM measurements in high- 



where Z is the partition function, |n), |m), |fc), En^m,k 
and Pn,ni,k are eigenfunctions and energy and momen- 
tum eigenvalues of the many-body Hamiltonian govern- 
ing the system.?/'''", ip are electron creation and annihi- 
lation operators. All operators are taken at the same 
spacial point 0. As follows from these expressions, the 



Tc superconductors, like in one may only conclude 
that the impurities are sufficiently far apart so that the 
Friedel oscillations from different impurities has room to 
develop without quenching each other. Therefore a rea- 
sonable starting point for the theory would be to study 
the Friedel oscillations of DOS from a single impurity. 
For a d-wave superconductor this was done in using 
BCS theory combined with a T-matrix approximation. 
The results are now widely used for the analysis of STM 
experiments on high-Tc superconductors giving very rea- 
sonable answers for the quasiparticle spectrum not far 
from the nodal points. Needless to say that since the cal- 
culation is based on a picture of well-defined weakly in- 
teracting quasiparticles, it cannot hold close to the antin- 
odal point where the STM measurements also yield very 
rich structure Q which requires a careful analysis. 

In general, Friedel oscillations of DOS are not very 
simple phenomenon and allow a straighforward interpre- 
tation only when there are well defined quasiparticles in 
the system. In the first order of perturbation theory in 
the impurity potential V one has 

(2) 
(3) 



relation between the Friedel oscillations and the excita- 
tion spectrum is not that straighforward. The situation 
is simplified only at small temperatures and when the 
many-body states in question can be approximated as 
quasiparticle ones. To see that let us consider ui > 0. At 
(3^^ = T = \n) = |0) is the ground state and Em — lo. 



5p{uj; r - ro) = ^ 5{uj + - Em){Anni + Bnm) 
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A.™ = E I ^^Pt^P-^^^ "^")] (fc|V-(0)|m)(m|^+(0)|n)(n|^(0)|fc)+c.c.l 
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Taking the Fourier transform of ^ in real space one fixes 
P/cm to be the external wave vector. Fixing Pkm rnay fix 
Ek and lead to singularity in the Friedel oscillation ampli- 
tude, but only in one-dimension and only if the eigcncn- 
ergies in question belong to quasiparticles. In D > 1 
equation Em — determines a surface in momentum 
space and the singularity is smeared by integration along 
this surface. In the same way the singularities are weak- 
ened when the eigenenergies En^k belong to continuum of 
states which cannot be parametrized by a single momen- 
tum. Therefore one is driven to the conclusion that in 
order to decipher the interference pattern concrete and 
model dependent expressions for DOS are necessary. 

In this paper I describe Friedel oscillations of local DOS 
for two strongly correlated systems - one-dimensional 
Mott insulator and one-dimensional Incommensurate 
Charge Density Wave (ICDW)/Superconductor. The 
motivation behind the calculations is twofold. First, 
these are strongly correlated systems which require non- 
perturbative approach. Nonperturbative results in this 
area predominantly concern Tomonaga-Luttiger liquids 
(see [4|). In the systems of choice where spectral gaps 
are generated dynamically and the long lange order is 
destroyed by gapless fluctuations the only available re- 
sult is [H]- Second, in my previous publication with 
Chubukov 0] we suggested that the dynamics of the 
antinodal regions in high-Tc superconductors is essen- 
tially one-dimensional. Therefore these calculations may 
be even relevant to the high-Tc problem after all. 

I will study the case when the spectral gap (be it 
the charge gap in Mott insulator or the spin gap in 
ICDW/Superconductor) is small compared to the band- 
width. The chemical potential is in the middle of the gap. 
The relative smallness of the gap enables me to employ 
field theory methods and bosonization. These methods 
together with the background information on the mod- 
els in question are described in may review articles and 
books, in particular in [3]. In the continuum limit both 
the Mott insulator and ICDW/SC are described by a 
universal Hamiltonian. The corresponding Hamiltonain 
density is a sum of two commuting parts governing dy- 
namics of the charge and the spin collective modes. For 
the Mott insulator we have 

n = rLc + rLs (4) 

—^icos{VSn(j)c) (5) 

= ^[K,id,e,f + Kr\d,<j)s)^] (6) 

where < 1, Kg and v^, Vg are the Luttinger parameters 
and the velocities in the charge and spin sector respec- 
tively and fi > is the coupling constant. Fields 4>aT0a 
obey the standard commutation relations: 

[9a{x),My)] = -^SabQH{x-y) (7) 



where Qh{x) is the Heaviside function. For ICDW/SC 
one has to interchange charge and spin indices. Then Kg 
becomes the Luttinger parameter in the charge channel. 
For 1/2 < Kg < 2 both 2kp charge susceptibility and 
pairing susceptibility are singular at zero frequency and 
temperature diverging as 

X{2kF,UJ^ 0) ~ Xpa^r{iO = 0) ^ T-^+V^^ 

(8) 

What order will eventually emerge when such one- 
dimensional systems are coupled together depends on 
whether Kg greater or smaller than 1 and to some ex- 
tent on the strength of the interchain interactions. Since 
both CDW and SC channels are singular, I do not distin- 
guish between the two and label the model ICDW/SC. 
However, to avoid confusion from now on I will discuss 
only the Mott insulator reserving ICDW/SC system for 
the final discussion. 

Model ([5]) is the sine-Gordon model; it is well studied 
and plenty of exact results are available, including results 
for its correlation functions. At Kc < 1 the spectrum has 
gap(s). The excitations always include solitons and anti- 
solitons carrying electric charge ±e (for ICDW it would 
be spin 1/2) and for < 1/2 also their neutral bound 
states (excitons). The asymptotics of the single electron 
Green's function is dominated by single soliton emission 
processes 

All spectral gaps are of the same order. The Mott gap 

is 

A(^/A)i/2(i-/<-c) 

where A is the ultraviolet cut-off and is probably of order 
of the bandwidth. 

For a clean Mott insulator DOS is featureless (apart 
from the spectral gap) and is given by ( see Eqs.(9,10) 
for = 2 in 0): 

p{u;)^{\u;\-mr'', 7 = ^( - Vv^)' (10) 

This is in contrast with band insulator; for Mott insulator 
with spin 1/2 electrons quantum fluctuations wipe out 
the singularity at the gap. We will see that impurities 
change this strongly enhancing the amplitude of Friedel 
oscillations at frequencies of order of m. 

Let us consider a single non-magnetic impurity coupled 
to the total density operator and placed at point xq . The 
most relevant contribution to the Hamiltonian density 
comes from the backward scattering and is described by 
the operator 

V = U{2kF) cos[\/2^0c(a;o)] cos[V2^<i>s{xn)] (H) 

where U{2kF) is proportional to the matrix element of 
the impurity potential. In what follows I will replace 
cos(-\/27r(/)c) by its vacuum expectation value. By doing 
so I neglect scattering of the massive solitons on fluc- 
tuations of the potential. It is justified when the bare 
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scattering potential is weak U{2kF) << A. The alter- 
native approach taken in Q was to replace the entire 
potential pT|) by a solid wall which is fine if the bare 
impurity scattering is close to the unitary limit and, in 
my opinion, substantially overestimates the contribution 
from such scattering when the potential is weak. Af- 
ter the above replacement the impurity scattering acts 
only in the spin sector where we have the boundary sine- 
Gordon problem: 



\5{x ~ xq)cos[V2tt4>s{xq)] (12) 

where A = C/(2/cf)(cos[V2^(?!)c]) - U {2kF){m/ K)^-'^ . 
This model is exactly solvable 11] and at the particu- 
lar value of the coupling Kg = 1 can be even reduced 
to the model of free fermions fl^. Another free fermion 
point (a trivial one) is Kg = 2. However, calculation 
of correlation functions of bosonic exponents is not an 
easy problem and in its entirety has not yet been solved, 
even at = 1. Some help comes from the fact that at 
Ks < 2 the impurity scattering potential scales to strong 
coupling and below the energy scale 

E* ^ A(A/A)2/(2-^=) 



can be replaced by Dirichlet condition 

4>s{xo) = 



(13) 



(14) 



This gives an easy way to calculate asymptotics of the 
correlation functions of the bosonic exponents using the 
method of images. In order to get a feeling for the 
magnitude of errors originating from deviations from the 
asymptotic regime, we will consider Friedel oscillations 
of the particle density 



{p{x) — po) = {co8{2kpx + V 27r(/)c) cos(v 27r(?!)s)) 
« cos(2fcFa;)m^'^/2(cos[\/2^(/)^(x)]) (15) 

where the charge cosine was replaced by its vacuum ex- 
pectation value. At i^s = 1 the average in the right hand 
side of (fTS]) was calculated in (see Eq. (3. 12)). Taking 
the Fourier transform of this formula we obtain 

{p{2kp + q)) ^ ^ ^ (16) 



E* 



One can check that a relative deviation of this function 
from its low-g asymptotic value exceeds 20 per- 

cent only at [q/E*) > 2.5. 

Now let us come back to our original task: calculation 
of oscillations of DOS. The bosonized expression for the 
fermionic operator is 



V.,(x) =e-i'=-^i?4x) 



Akpx 
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(17) 



(18) 



where rja- are the Major ana zero modes (Klein factors), 
ao is the lattice constant and cr = ±1. 

The single-particle density of states is related to the 
single-particle Green's function at coinciding spacial 
points: 



G(r; X, x) — Ggmooth + GoscU 

Gsmooth = {{R{t,x)R+{0,x))) 

{{L{t,x)L+{0,x))) 



G. 



oscil 



,21fei 



"((L+(t, x)R{0,x)))+H.c. 



(19) 

(20) 
(21) 



where r is Matsubara time. Let the reader note that 
Goscii = in a clean sample. 

The operators R, L and their Hermitian conjugate fac- 
torize into charge and spin parts (|18p governed by dif- 
ferent Hamiltonians ([U [T2|) respectively. Since the im- 
purity contributes mostly to the spin Hamiltonian, for 
the charge sector one can use the results obtained in 
Q. Then for the oscillatory oscillatory part (pij) of the 
Green's function we obtain 



Goscii = am^^^Ko{mT)S{T, x - xq) 
S{t,x- xo) = 



(22) 
(23) 



where a is a numerical coefficient. The Kq function in 
comes from the charge sector (the corresponding cor- 
relator was calculated in 0]). Therefore the fact that in 
the clean system object (dJ) vanishes altogether is due to 
the spin sector remaining critical. The latter means that 
its right and left sectors do not couple to each other and 
S = 0. The scattering from the impurity(ies) connects 
the right and left sectors and S{t, x) ((23)) is no longer 
zero. Thus impurities reveal the parity violation which 
exists in the Mott insulator already in the clean case, but 
is not directly observable for a lack of a suitable local op- 
erator. 

Function ([23| can be represented as 



S{t,x) 



1 VgT ^ 

jT( ,xE /Vs) 



(24) 



where d ~ {K+\/K)/^ and T{y, z) is a scaling function. 
In the strong coupling limit (that is at the distances or 
times much larger than l/[£^*]) the impurity is substi- 
tuted by the condition ([14]) . This condition is equivalent 
to the conditions 

cj)s{x-xo) ^ -'(j>s{xo-x), 9s{x-xo) = 9{xo-x) (25) 

Then correlation functions of bosonic exponents are cal- 
culated by the method of images, like in electrostatics. 
As a result in the strong coupling limit we obtain univer- 
sal (that is independent of the bare impurity potential) 
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asymptotics: 



frequency. This is due to the fact that the electron spec- 
trum is hnearized and the chemical potential is exactly in 



Goscii = acos[2kF{x - xq)] x (27) 

mi/"i^o(mr)|r|(^--i/^-)/4 
|(:r - xo)/vs\^^/i[T^ + iix - xo)yv^,]^^/» 

where factor Z was calculated in [§] and a 1. We 
perform a double Fourier transformation of (I27p , both in 
space and Matsubara time and then perform the neces- 
sary analytic continuation iw — > w -I- iO. The answer can 
be written as 



dp{uj > 0; 2kF + q) 

r2{uj — m) /vs\q\ 

-l+Ks/2 



(28) 

dyA{i,-y\q\vj2)f{y) 

f{y)= I dz(y2-,2)-i+AV8|i + ,|-i+K./4 (29) 

-y 



Here 



-1/2-0 



(lu / m) 



a/2 



F (-a/2, 1/2 - a/2; 1/2 - a; 1 - (m/Luf) (30) 

where a — {Kg — 1 /Kg) /4 and F{a, b; c; x) is the hyperge- 
ometric function. Numerical evaluation of these integrals 
show that that at 1/2 < < 2 the Fricdel oscillations 
do not disperse and can be approximated by 

5p{LU, q + 2kp) - \q\-^+^^/^g{uj/m - 1; K^) (31) 

where function Q{x; Kg) is depicted on Figs. 1,2 for Kg = 
1,2. 




FIG. 1: Function g{x-l). 

Thus we have found that impurities in ID Mott insu- 
lator (or ICDW/SC) produce 2kF Friedel oscillations of 
DOS. These oscillations occur at frequencies larger than 
the Mott gap (the soliton gap) and at T = their am- 
plitude has a maximum ed to ~ 3m and a strong singu- 
larity at 2kF- Near the singularity the Fourier transform 
of the DOS is essentially dispersionless. In our calcula- 
tions DOS (including the Friedel oscillations) is even in 




FIG. 2: Function e(x;2). 



the middle of the gap. The latter feature can be trivially 
corrected by a chemical potential shift. 

For Mott insulator with SU(2) symmetry in the spin 
sector Kg — 1 and a — 0. For ICDW/SC spin and charge 
sectors s and c are interchanged; then Kg corresponds 
to the Luttinger parameter in the charge sector and its 
value is not fixed by any symmetry. The phenomenology 
developed for the cuprates in @ suggests Kg 2. Note 
that at Kg — 2 operator pT|) becomes marginal and E* 
increases exponentially when the value of Kg approaches 
2 from below. 
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